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SUMMARY 


Hansen’s method of partial anomalies is described and 
his application of the method to Encke T s comet is discussed. 
The principle innovation of the method is to divide the el- 
lipse of the perturbed body into a small number of seg- 
ments. In each segment an independent variable closely 
related to the conventional anomalies is selected in order 
that the series representing the perturbations is strongly 
convergent within the segment but is invalid outside of it. 
A different set of series is obtained for each segment of the 
ellipse. In the case of high eccentricity, these series can be 
made much more convergent than the conventional series. 

This method may be applied to the determination of the 
perturbations of artificial satellites and space probes with 
highly eccentric orbits. The’TMP and AIMP fT type satellites 
and probes to Venus, Mars, or Jupiter could be handled by 
this theory. In fact, any celestial object, natural or arti- 
ficial, with an eccentricity greater than about 0.50 may be 
very efficiently dealt with by means of the method of partial 
anomalies. This method has the great advantage of being 
rapidly convergent for all eccentricities and of being 
semianalytie. 
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AN EXPOSITION ON 

HANSEN'S METHOD OF PARTIAL ANOMALIES 

by 

George E. McCluskey, Jr. 

Goddard Space Flight Center 


INTRODUCTION 

The perturbations which one planet produces on another or on a comet depend on the mean 
anomalies of the two bodies. The perturbations may be expressed in the form: 



where p. and p are functions of the elliptic orbital elements of the two bodies and of i and i ' 
which are integers, i extending from 1 to •> and i' from - - to + The mean anomaly of the per- 
turbed body is g while that of the perturbing body is g ' . 

If one of the bodies has an orbit of high eccentricity and is highly inclined with respect to the 
other, the above series will converge very slowly, if at all. In the case of a comet the radius vector 
may at one time be much smaller than that of the perturbing planet and at another time much 
larger. Thus, a series using one of the conventional anomalies as its argument may converge in 
one part of the orbit but will fail to do so in another part. 

In an effort to correct this situation, Hansen (Reference 1) has chosen a new variable to rep- 
resent the motion of the perturbed body in place of the mean anomaly. This new variable is called 
the partial anomaly. The partial anomaly, when applied properly, will lead to very rapidly con- 
verging series with respect to the perturbed body for any value of the eccentricity less than unity. 
The method may even be applied to parabolic and hyberbolic motion. 

Hansen introduces two partial anomalies, one known as the inferior anomaly and the other as 
the superior anomaly, and also an intermediate anomaly. By means of these anomalies the ellipse 
is divided into two, three, or any number of segments. Each segment is represented by one par- 
ticular partial anomaly. The equations of the ellipse and the equations of elliptic motion are ex- 
pressed in terms of the partial anomalies, as are the perturbations. Each segment has its own set 
of series for the perturbations. A given set of series cannot be used outside of its particular seg- 
ment but will converge rapidly within its segment. 
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I. THE PARTIAL ANOMALIES 


L.et: 

f = True anomaly 
u = Eccentric anomaly 
k = Inferior partial anomaly 
kj = Superior partial anomaly 
a = Semimajor axis 
c = Eccentricity 
r = Radius vector 

n = Mean motion (Kepler’s third law) 
t = Time 

Define k and k : by: 


l 

sin 2 u - t ( cos Xsink + sinX) 


cos 2 f tj (cos Xj sinkj -sinXj) (i) 


where: 


r r ' + r" - 2a( 1 -e)] 1/2 

t: “ L 4ac 

[~a(r'+r")(l + e)~2r'r" 1 

- * I' 4er' r - (1 “ e) J 

f r" ~a(l - 

Lr' -a(l - 


1/2 


tan (45 - X) - 


■_o 

■e). 


tan 


(4S _ x ) , R l+e >- . r .: i: 

v 11 |_a(l+e)-r" r'_ 


1/2 


r' and r" being any two radius vectors of the ellipse. 
Equations (2) show that: 


(2) 


0 < t < 1 

0 < c, < 1 


2 


The equations of the ellipse are 


7 7 

"4 


7 T 

7 


< X < +^r 


r ' - a( 1 - e cos u ' ) 
r " - a( 1 - e cos u" ) 

or: 

r . = a ( 1 ~ e 2 ) 

1 + e cos f ' 

r „ = a (i - e2 ) 

1 + e cos f " 

Using Equations (3) and (4), we may write Equations (2) in the following forms: 


or: 


„ . 2 1 ; 

u + s in o - u 


, 1/2 



cos 2 ~2 f w + cos 2 2 1 

e l ■ 

2 


1 


sin 9 u" 

tan (45 -X) = x 


sin ^ u 


1 r, 


COS -77 f 

tan (45 

-*.) = - 1 f „ 


cos 2 f 


f cos X 


1 

sin 2 u" 


+ sin 2 u ' 

2 


e sinX 


1 

-sin 2 u" 

2 


( 3 ) 


(4) 


(5) 


3 


€ j COS Xj 


cos 2 f 1 ” cos f " 


COS f ' + < 


€j sin Xj 


sin 2X 


. 2 1 , . j 1 „ 

sin z j u ” sin ^ u 


(r 1 ~ r" ) (l - e 2 ) i 
4e r 1 r " e t 2 


cos 2 j f“ " cos 2 -j f 1 


The equations of elliptic motion are: 


Since, by Equation (1): 


r - a( 1 - e cos u) 


r cos f - a(cos u - e ) 


rsinf - avl - e 2 sinu 


n dt - - du 


sin o u e 


[cos X sin k + s in X] 


we write Equations (8) in terms of k. 


r - a(l-*e+et. 2 +et 2 sin 2 x) + 2aet 2 sin2Xsink - act 2 cos 2 X cos 2k 


cos f - a(l-e-e 2 -e 2 s in 2 X ) - 2at 2 sin 2X sin k + a< 2 cos 2 X cos 2k 


rsinf - 2at y 1 - e 2 (cos X sin k + sin X) A 


r ( 2t cos X cos k dk ' 


A ( 1 - <- 2 sin 2 X - t. 2 sin 2X sin k - * 2 cos 2 X s in 2 k) 1 
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On the other hand, the elliptic equations may be written as: 


_ a(l - e 2 ) 

1 1 + e cos f 


n dt 



and by Equations (1): 


( 10 ) 


cos 2 f - (cosXj sinkj -sinXj) 

we write Equations (8) in terms of k x . 


(l - e - ee j 2 + ee 2 sin 2 Xj) - 2ee j 2 sin 2Xj s in kj - ee 2 cos 2 Xj cos 2k j 

a(l- e 2 ) 


cos f 

(e j 2 + Cj 2 sin 2 Xj - 

1 ) - 2e j 2 sin 2Xj 

sinkj - ej 2 cos 2 Xj cos 2kj 

sin f n 

2e j (cos Xj sin k j 

- sin Xj) B 


n dt - 

r 2 

cos Xj cos k j dkj 


a 2 |/l - e 2 

B 


B = 

( 1 - e j 2 sin 2 Xj + e 

2 sin 2Xj sinkj - 

6 j 2 cos 2 X j sin 2 k j ) 1/2 


(ID 


The equations of elliptic motion are now expressed in terms of the partial anomalies. When a 
and e are known, we can calculate e, e l9 X, and X l for any pair of values of r ' and r". 

The first of Equations (9) gives: 


r - a(l-e+ee 2 +ee 2 sin 2 x) + 2aee 2 sin2Xsink - aee 2 cos 2 X cos 2k 


Consequently: 


dr 

dk 


2aee 2 (sin 2X cos k + cos 2 X sin 2k) 


The maxima of r are at: 


dr 

dk 


0 
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which gives: 


k = 90° 


and 


k = 270° 


These values give: 


r = r' at k * 90° 

r = r" at k = 270° 


The minimum value of r is at: 


sink - tanX 


which gives: 


r - a( 1 ~ e ) 


Similarly, using the first of Equations (11): 


kj = 90° at r “ r ' 

kj - 270° at r =■ r" 

where r' and r" are now minima of r. The maximum of r is at: 


r - a( 1 + e) 

The ellipse has been broken into two parts. The partial anomaly k represents the portion of 
the ellipse from perihelion to r - r' on one side of the major axis and from perihelion to r - r" on 
the other side. The partial anomaly represents the ellipse from aphelion to r = r' on one side 
of the major axis and from aphelion to r = r" on the other side. At r - r", k = 270°. As r de- 
creases, k increases. At perihelion, k = 0°. For 270°5k < 360°, we have the true anomaly, f, 
greater than 180 . After passing perihelion, r increases as does k until at r r' we have k = 90°. 
For 0 k 90 , we have f 180°. If we allow k to increase beyond 90° we retrace the same portion 
of the ellipse in the reverse order. 
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At r - r # we have k x = 90°. As r increases, so does k x . At aphelion we have k x = 180° and 
finally at r = r" , k 1 = 270°. If k x increases beyond 270°, the same portion of the ellipse is retraced 
in the reverse order. 

We see that k, the inferior anomaly, represents the portion of the ellipse containing perihelion 
and cannot represent the portion containing aphelion unless r f = r" = a(i +e) , in which case we 
have not really divided the ellipse at all and k is equivalent to the eccentric anomaly. Similarly, 
k x represents the portion of the ellipse containing aphelion. 

By selecting one of the points of separation, r' or r", at perihelion or aphelion, we may divide 
the ellipse into three or four segments. Consider how we divide the ellipse into four parts. 

Let: 


r " - a( 1 - e) 

k - 2Kj - 90 

E = Y2 e 

Equations (9) become: 


r 

= a (l 

3 \ 1 

- e + ^ eE 2 j - aeE 2 cos 2K X + ^ aeE 2 cos 4K 

r cos f 

= a (l 

~ e + 4 E 2 ] + aE 2 cos 2Kj - -4 aE 2 cos 4K r 

r sin f 

~ aE j 

/l-e 2 (l-cos2K 1 )y'l-E 2 sin‘ , K 1 

n dt 

r 

2E sin 2K, dK, 

*a 

1^1 - E 2 sin 4 K. 


Equations (12) show that at perihelion K x = 0°. As K x increases, r increases. At K x = 90°, 
we have r = r ' . Thus, these equations represent the portion of the ellipse from perihelion to r = r ' . 
The quadrant is defined by noting that Equations (12) give f < 180°. For values of greater than 
90° the same segment of the ellipse is retraced. 

Let: 

r " = a ( 1 + e ) 

kj = 2K 2 - 90° 

El = 


7 


Equations (11) become: 


l (l - e + 4 - eEj 2 ) - eEj 2 cos 2K 2 + \ eE 2 cos 4K 2 

7 = a(Y-~?) 

cos f = ( 4 “ Ej 2 - l) - E 2 cos 2K 2 + -4 Ej 2 cos 4K 2 

sin f = Ej ( 1 ~ cos 2K 2 ) y^l - Ej 2 s in 4 K 2 

r 2 2Ej sin 2K 2 dK 2 

n dt = • • ■■ — 

a 2 1 - e 2 ^1 - Ej 2 sin 4 K 2 

Equations (13) show that at r = r' , K 2 = 90°. As K 2 increases, so does r until at aphelion, 
r = a(l +e), we have K 2 = 180°. For values of K 2 outside of the range 90°£K 2 < 180°, the same 
segment of the ellipse is represented. We also see that f £ 180°. 

Let: 


t ' - a(l + e) 
k x = 2K 3 + 90° 

Ei = a ej 

Equations (11) become: 

1 (l - e + 4 : eE t 2 ) - eE t 2 cos 2K 3 + 4 eEj 2 cos 4K 3 
7 = a (Y ~ e 2 ) 

cos f = (| Ej 2 - 1 ) - El 2 cos 2K 3 + ^ E , 2 cos 4K 3 


sin f = Ej (cos 2K 3 - l) ^1 - E , 2 sin 4 K 3 


r 2 2Ej sin 2K 3 dK 3 

n dt y • — 

a 2 ^l - o 2 ^1 - Ej 2 sin 4 K~ 


(14) 


Equations (14) show that at aphelion, K 3 = 180°. As K 3 increases, r decreases. At k 3 = 270°, 
we have r r" . We also have f > 180°. 
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Finally; 


r' = a(l-e) 
k - 2K 4 + 90° 

E = € a 


Equations (9) become: 


r - a (l - e + ^ eE 2 ) - aeE 2 cos 2K 4 + -4 aeE 2 cos 4K 4 
r cos f = a(l-e“ 4 ;E 2 ]+ aE 2 cos 2K 4 - -4 aE 2 cos 4K 4 


r sin f - aE ^1 - e 2 ( cos 2K 4 r 1 ) ^ 1 " E 2 sin 4 K 4 


n dt 


r 2 E sin 2K 4 dK 4 


^1- E 2 sin 4 K 4 


(15) 


When K 4 = 270°, r = r" and as K 4 increases r decreases until at perihelion we have k 4 = 360° 
Also, f > 180°. 

Table I summarizes the results. 


Table I 


Partial Anomaly 

Range 

Range of r 


0° - 90° 

perihelion to r ' 

K 2 

90° - 180° 

r ' to aphelion 

K 3 

180° - 270° 

aphelion to r " 

K 4 

270° - 360° 

r " to perihelion 


The ellipse has been divided into four parts with the points of separation at r D a(i-e), r = r', 
r = a(l + e) , and r = r" . 

Suppose we omit the point of separation at perihelion and take the three points of separation 
at r = r 1 , r = a ( 1 +e), and r = r" . Then Equations (9) using the inferior partial anomaly k 
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represent the portion of the ellipse from r" to r' containing perihelion. We have: 

k = 270° at r = r" 

k = 0° at r - a(l-e) 

k = 90° at r = r' 

The portion of the ellipse from r' to r = a(i + e) is represented by Equations (13) using K 2 . 
We have: 


K 2 = 90° at r = r’ 

K 2 = 180° at r - a( 1 + e) 

The portion of the ellipse from a(l + e) to r" is represented by Equations (14) using K 3 . We 
have: 


K 3 = 180° at r - a( 1 + e ) 

K 3 = 270° at r = r" . 

Thus, the ellipse has been divided into three sections. 

The partial anomalies k and k t and the related anomalies K , , k 2 , K 3 , and K 4 require that the 
points of separation be on opposite sides of the major axis or at the end points of the major axis. 
We cannot take two points of separation on the same side of the major axis. Since it may be nec- 
essary to select two or more points of separation on the same side of the major axis, we define the 
intermediate anomalies which permit the number and position of the points of separation to be 
completely arbitrary. 

Let r, and r 2 be any two radius vectors of the ellipse on the same side of the major axis. For 
definiteness, letr 2 > r r Define the intermediate anomalies L, and L 4 by: 


and define E by: 


cos Kj - i cos Lj 


cos K 4 - £ cos L 4 

(16) 

rr - a( 1 - e >1 1/2 


[ 2ao J 

(17) 
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The modulus l will be defined later. Equations (12) become: 


r - a \l -e + eE 2 (2- 2l 2 + f ^ 4 )J - aeE 2 { 2 I 2 - V) cos 2L X + \ aeE 2 V cos 4Lj 
r cos f = a jl - e - E 2 (2 - 2^ 2 + 4; *i 4 )j + aE 2 ( 2^ 2 - ^ 4 ) - 4: aE 2 V cos 4L X 

r sin f = B' aE^f-e 2 [(2 -l 2 )- 1 2 cos 2 lJ 

r 2&£ 2 sin 2Lj dL l 

ndt = - • — 

a 3' 

where: 


(18) 


B # 



<e 4 ) 


cos 2 Lj 


1 1 1/2 
g E 2 't 4 cos 4Lj 


Equations (15) become: 


r 


r cos f 


a|l - e + eE 2 (2 - 2 ^ 2 + 4 ^ 4 )j ~ aeE 2 ( 2 ^ 2 - ^ 4 ) cos 2 L t + 4 - aeE 2 Z 4 cos 4L 4 
a j^l - e - E 2 (2- 2-t 2 + 4 ^ 4 )] + a£2 ( 2 ^ 2 ~ V) cos 2L 4 - \ aE 2 V cos 4L 4 


r sin f 


aE |/l - e 2 i 2 cos 2L 4 - (2 - 1 2 )] C' 


n dt 


r 2Ef , 2 sin2L 4 dL 4 
a c « 


(19) 


where: 


C' - [l-E 2 (l--C 2 + |r) +E 2 (-€ 2 - ^ cos2L 4 - ^E 2 r cos4L 2 ] 1/2 

The expressions for Lj are used if the values of f t and f 2 , corresponding to r t and r 2 , are 
less than 180°. If f x and f 2 are greater than 180°, the expressions for L 4 are used. Equations (18) 
give r - r 2 atLj = 90° while Equations (19) give r = r 2 at l 4 = 270°. Since r 2 is the maximum 
value of r to which Equations (18) and (19) apply, other values of l must pertain to values of r less 
than r 2 . Since the equation for r sin f in (18) and (19) shows that sin f cannot change sign, all of 
the values of r represented must lie on the same side of the major axis whatever l may be. 
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Obviously, we wish r = r, to be the minimum value of r represented. We thus define the 
modulus l of Equations (16) such that at r = r, , we have L = 0°; that is we have l, = 0° or l 4 = 360 
If L = 0° the first of Equations (18) or (19) gives: 


This gives: 


a[l - e + eE 2 (2 - U 2 + 2-t 4 )] 


'C 2 


r, -a(l -e)-j 1/2 
r 2 ~ n ( 1 " c )_ 


( 20 ) 


Define Cl) by: 


'r, -a(l -e)-j 1/2 
tan ( 45” - co) * [^-^1 - e ) 


( 21 ) 


Then: 


-C 2 


sin n- j- 

cos (45° _ o ) 


( 22 ) 


This definition of l gives r r x for L x = 0° or l 4 = 360°. The value r = r 2 occurs for L x = 90 
or L 4 = 270°. Values of L x outside of the range 0°< L l < 90° or values of L 4 outside of the range 
270° ^L 4 - 360° simply retrace the same segment of the ellipse between r - r t and r r 2 . Con- 
sequently, L x represents the segment of the ellipse from r r r x to r - r 2 if f < 180° and l 4 rep- 
resents the segment of the ellipse from r^r 1 to r =r 2 iff> 180°. Both r x and r 2 must be on the 
same side of the major axis if we are to apply the intermediate anomalies. If they are not, of 
course, we would apply the partial anomalies k and k , . 

We may also define l 2 and l 3 . 


COS K 2 J COS L 2 

COS K 3 ~ 'i x cos L 3 (23) 


Define: 


a( 1 4 e ) - r J 

2 c r j (>-<•) 


( 24 ) 
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We could then write Equations (13) and (14) in terms of l 2 and l 3 , respectively. We then 
define ^ as: 

sin a) j 

^1 " cos (45° ^ (25) 

where: 


tan (45° - oj x ) 


a(l + e) - r 2 
a(l + e) - r , 



(26) 


We would then find that r = r x for L 2 = 90° or L 3 = 270°. The value r = r 2 occurs for h 2 = 180° 
or l 3 = 180°. Values of L 2 outside of the range 90°<l 2 < 180° or values of L 3 outside of the range 
180° < L 3 < 270° will retrace the same segment of the ellipse between r = r x and r = r 2 . Conse- 
quently, L 2 represents the segment of the ellipse from r = r x to r = r 2 if f < 180° and L 3 represents 
this segment if f > 180°. 

As an example, consider six points of separation at r q a(i - e) , r = r x , r = r 2 , r = a(l + e), 
r = r/ , and r - r 2 • Figure I shows the situation. The application of the various partial anomalies 
is given in Table II. 


We could divide the ellipse into more than six parts by defining more intermediate anomalies 
and moduli in the same manner. 



L 4 

Figure I— Division of the ellipse into six sections. 
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Table II 


Partial Anomaly 

Range 

Range of r 

K, 

0° — 90° 

perihelion to r t 

L , 

0° — 90° 

r i t0 r 2 

k 2 

90° — 180° 

r 2 to aphelion 

K 3 

180° — 270° 

aphelion to r 2 ' 


270° — 360° 

r j to r / 

k 4 

270° - 360° 

r / to perihelion 


In applying the partial anomalies to the calculation of perturbations, the various radicals 
appearing must be developed in a Fourier series of multiples of the partial anomaly. This devel- 
opement will be made so that very convergent series for the perturbations are obtained. The 
manner in which we select the number and position of the points of separation will depend on the 
specific problem. They will always be selected to make the various moduli, t, c lf E, E t , l 9 or l v 
as small as possible since the smaller the modulus the more convergent the series obtained will be. 

Consequently, by dividing the ellipse properly, the series representing the perturbations may 
be made to converge rapidly. A different series is used in each segment of the ellipse. If the el- 
lipse is divided into n parts, we will have n sets of series for the perturbations. But the convergence 
will be much more rapid than that of the series which would represent the whole ellipse using the 
mean anomaly or one of the other conventional anomalies. 
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II. DEVELOPMENT OF THE RADICALS 


In Equations (9), (11), (12), and (18), radicals of the following form appear: 


(l - e 2 sin 2 k) ±1/2 

(27a) 

(l -E 2 sin 4 k)* 1/2 

(27b) 

(l + sin 2 K ) i 1/2 

(27c) 

(A + Bsink + Csin 2 k) ±,/2 

(27d) 

(A' +B' cos 2L + C' cos 4L) 11/2 

(27e) 


We shall be able to write all of these in a form involving only factors of the form of Equation 
(27a). Consequently, consider the expansion of Equation (27a). Let: 

(l - e 2 sin 2 k)“ 1/2 = a Q - 2ol 2 cos 2k + 2a 4 cos 4k - *•* (28) 

where: 


a 


2i 



cos 2ik dk 
(l - e 2 sin 2 k) 1/2 


(29) 


The plus sign is for i even and the minus sign for i odd. 

The following recursion formula may be obtained by integrating by parts or by developing the 
integrand as a power series in k and comparing terms for three successive values of i. 

0 - (21 -1) e 2 a 2 ._ 2 - 4i(2- 6 2 ) a 2i + (2i + 1) e 2 a 2i+2 (30) 

If: 

(l - e 2 sin 2 k) + 1/2 = ~ 2/3 2 cos 2k + cos 4k - (31) 

then: 


P 


2 i 


8i ( a 2i-2 a 2i+2) 


(32) 
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except: 


- (* " 2 fi2 ) a o - \ a 2 ( 32 ) 

Equation (32) may be obtained by recognizing that: 

(l-e 2 sin 2 k)~ 1/2 (l-e 2 sin 2 k) + 1/2 - 1 

The calculation of a 4 , a 6 , etc. from a 0 and a 2 by means of Equation (30) is inaccurate as er- 
rors in the last significant figures of a 0 and n 2 cause larger and larger errors in a 2i as i in- 
creases. To circumvent this problem, the following scheme is provided. 

Define: 


Equation (30) becomes: 


4i (2 - e 2 ) 

(2i - l) ! 2 P21 

(2i-l) 2 / ^ V 

I6i(i - 1) ^2 - e 2 / 


0 


y 2i r 


2 i + 2 2 i + 2 


(33) 


(34) 


Equation (33) show that as i increases, the quantity y 2i+2 -y 2i decreases. At i - 

^2i 0 ^ 2 i + 2 ^oo a ^ S0: 



Equation (34) becomes: 
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or: 


y - = - 2 k £ / - [(2-€ 2 )±2 yr^T 2 ] (35) 

The (-) sign in front of the radical must be used as the (+) sign leads to divergent results. 

Define: 


6 = s in \p 


Equations (2) lead to: 

( i ' 
cos \p - tan (45 - j X 


Equation (35) becomes: 

2 1 

y a - sec 2 2 x 


Equation (38) is exact for i - The definition of 'P allows us to write Equation (33) as: 


2i - 1 

P 2i - — 4T sm X 7 2i 

2i - 1 2 „ 

^21 - 16i( i - 1) Sln X 


(36) 


(37) 


(38) 


(39) 


We proceed as follows. Select some large value of i = i f . Compute k 2i for i = 2, 3, 4, ... i f 
by Equations (39). Then, assume y 2i = y a . We calculate y a by Equation (38). Then, Equation (34) 
is written as: 


■X2i-2 1-X 2i y 2i 

For i = if, use y 2i = y ^ and calculate r 2i _ 2 . Then calculate y 2i _ 4 , y 2i _ 6 , • • • , y 2 . The first of 
Equations (39) is used to calculate p 2 , p 4 , • ■ • p 2i . By Equation (33): 

a 2i = a oP 2 P4’"P2i (40) 

In this manner we compute all of the a's from a 0 . The equation y 2i = y w is only approximate 
for finite i. The accuracy of a given calculation can be estimated by recalculating the >'s for a 
larger value of i f and noting the effect on the values of the a's. 
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We may improve the approximation to y 2i in the following manner. Recall that: 


0 


r 2i y 


2i '2i + 2 v 2i + 2 


and: 

(2i - l) 2 

X 2i 16 i ( i " 1) Sln X 
or: 

(2i + l) 2 

X 2i+2 “ 16i(i + 1) s,n X 

Using this expression for ^ 2i+2 and sec2 1/2 X for y 2i+2 , we have: 



- l 

r (2i + l) 2 , 1 


Vn - 1 ■ 4i(i + 1) s,n 2 X 


L J 



Developing this expression and dropping powers of sin 1/2 X higher than the second, we have: 


y 7i - sec* 


2 X j^ 1 + 4i(i + 1) sin2 2 X j 


(41) 


Since X< 45°, the second term is always less than 1/1000 of the first term for i = 10. Conse- 
quently, this approximation is sufficient for all cases of interest for moderate values of i. 

It remains to calculate a 0 . Equation (29) gives: 


a 


o 


z * 77 /2 

2 _ f dk 

77 J 0 (l-c 2 sin 2 k) 1/2 


or: 


a 


o 



(42) 


f (e, v/2) is the complete elliptic integral of the lirst kind. It has been extensively tabulated. It 
may be calculated from the following series if e 2 is sufficiently small. 


9 



<» r 

V" ( 2n ) ! 

Z_> 2 2 " (n! ) 2 


,, 2n 


(43) 
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This series is easily obtained by expanding (l - e 2 sin 2 k)“ 1/2 by the binominal theorem and in- 
tegrating the series term by term. The series given by Equation (43) will converge fairly rapidly 
if e 2 < 1/2. 

Consider now the form of the radical given by Equation (27b). We may write: 


and: 


(l - E 2 sin 4 k)~ 1/2 = (l + E 2 sin 2 k)" 1/2 (l - E 2 sin 2 k) 1/2 

" + 2V 2 cos 2K + 2A. 4 cos 4K + ••* 


(44) 


(1 + Esin 2 K)~ 1/2 = e 0 + 20 2 cos 2K + 25 4 cos4K + ••• (45a) 

(l - E sin 2 k)' 1/2 = r) 0 - 2t] 2 cos 2K + 2t} 4 cos 4K - ••• (45b) 

By multiplication of series: 

\ = ^0 9 0 - 2t ^2 6 2 + 2t ?4 9 A ~ ' " 

“S = ^2 9 0 _ N + ^ 4 ) 6 2 + (v 2 + ’O 0 4 _ ‘ ‘ 

X 4 = ^4^0 ‘ ( T '2 +T ? 6 ) 0 2 + (V> + V>) 0 4 " (46) 


Equation (45b) is in the same form as Equation (28) and may be handled in exactly the same 
manner if we put: 


yE = sin \p 


^ 2 i - ^0 P 2 P 4 ‘ ' ' P2i 


(47) 


We also have: 


so that: 


^0 



dK 

(l - E sin 2 k) 1/2 



(48) 


Thus, we may calculate the Vs. 
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Now consider Equation (45a), Define: 


K - 90° - P 

Then: 


l + Esin 2 K - ( 1 + E) ^1 - ^ + £ sin 2 

(49) 

This is now of the same form as Equation (28) and is dealt with accordingly, 
late the £ f s and by Equations (46), the \ f s. 

Equation (27c) may be changed to the form: 

We may now calcu- 

(l + sin 2 K)" 1/2 = + oj 2 cos 2K + , 4 cos 4K + 

Define: 

K 90° - P 

Then: 

(50) 

(l + sin 2 k)“ 1/2 - (l ^ cos 2 P)" 1/2 = 2^1 - ^ s i n 2 P^ 

(51) 

This is again of the same form as Equation (28) and we may calculate the u T s. 
Consider the form of Equation (27d). Write: 

(A + B sink -C sin 2 k) 12 - + 2\jSink + 2k 2 cos 2k + ** 

(52a) 

(A + B sink - C sin 2 k) + 1/2 - + 2/± l sin k + 2/^ 2 cos 2k + ** 

(52b) 


we will have: 




i 77 


i r _ 

" J. («- 


cos ik dk 


B s in k - C sin 2 k) 1 : 


r n 

“77 (A + B sink - C sin 2 k) 1/2 cos ik dk 

Jo 


(53) 
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for i even. For i odd: 


Since: 


we have: 






1 

7 7 


1 . 

77 




sin ik dk 

+ B sin k - C sin 2 k ) 1/2 


+ B sin k - C sin 2 k) 1/2 


sin ik dk 


(A +B sink -C sin 2 k) 1/2 (A + B sin k - C sin 2 k) + 1/2 = 1 


(54) 


4i (^i-2 ^i + 2 ) 4i (^i-i + ^i + i) 

and: 

^ = ^(^i- 2 '^i + 2 ) + l"(^i-i + ^i + i) (55a) 

for i even and odd, respectively. The exception is: 

^0 = (a - ^c)\ 0 + BX, + c\ 2 (55b) 

The coefficients A, B, and C will be given by: 

A = l~e 2 sin 2 X 
B = e 2 sin 2X 

C = e 2 cos 2 X (56) 

Then: 

(A + B sin k - C sin 2 k)~ 1/2 = (l - £ 2 sin 2 X + e 2 sin 2X sin k - e 2 cos 2 X sin 2 k) " 1/2 

- (l - « sin X + € cos X sin k) ‘ 1/2 (l t e sin X - e cos X sin k) ' 1/2 (57) 
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Define: 


(l - e sin X + e cos X sin k) 1/2 - 6 0 - 2^ sink 

and: 

(l + e sin X - € cos X sin k)" 1/2 = r) Q + 2 77 t sin k 

Equation (52a) and Equations (58) give: 


^0 = 

V 0 ^0 

~ ™ > 


+ 2 e 2 

v 2 ~ 

2^3 

^3 + 

- 

Vi 0 O 

- 0 X 1 

fro 

+ r ^ 2 ) 

+ *2 

hi + 


^2 = 

v 2 0 o 

- 0, 


+ ^) 

+ 0 2 

(^0 + 


-*3 = 

v 3 0 o 

- *i 

(^2 

+ v A ) 

+ 0 2 

hi + 

^ 5 ) 

\ = 

r >*°o 

- *i 1 

(^3 

+ r >s) 

+ 


' T '«) 


To determine 0. let: 

k 


and to determine let: 


k 


Then: 


1 - t sin X + e cos X sin k 


and 


1 t e sinX " ecosXsink 


26 2 cos 2 k + 26 z sin 3 k + 2 0 4 cos 4 k - * * * (58a) 


2 r} 2 cos 2 k - 2t) 3 sin 3 k + 2 t ? 4 cos 4 k + * * * (58b) 


20,v A - ■■■ 


- 

(^2 +7 U) 

+ e 4 


+ 77 5 ) - ••• 


- *3 

K +7 5s) 

+ ^4 

hi 

+ r , 6 ) - ... 


" 0 3 

ho ■ 

h7 J«) 

+ ^4 

hi 

+ v 7 ) - ••• 


" *3 

hi 

^ 7 ) 

+ o A 

ho 

+ 7?8 ) _ ... 

(59) 


21 + 90° (60a) 


2V - 90° (60b) 


(l - £ sin X • £ cos x) (l - f 2 sin 2 -t!) (61a) 


Jl t £ sin X - i cos x) (l - f ' 2 sin 2 f. ') (61b) 
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where: 


f 


f ' 


2 e cos X 


1 - e sin X + e cos X 


2e cos X 


1 + e sin X - e cos X 


1/2 


1/2 


Thus, we have again obtained the form of Equation (28) and may compute 0. and 77 . . 
follow from Equation (59). 

Equations (58) give: 


/ 2 

2 f dk 

77 Jo [( 1 - € sinX) + e cos X sin k] 1/2 

2 C dk 

77 J [( 1 + e sin X) - e cos X sin kj 1/2 


Define: 


a = 1 - e sin X 


b = e cos X 
a ' = 1 + e sin X 


Since -45° < X< +45°, a, b, and b' are always greater than zero. Then: 


^0 



dc 

(a + b sin k ) 1/2 


dk 

(a' - b sin k ) 1/2 


Equation 288.00 of Byrd and Friedman (Reference 2) gives for a >b > 0: 


f 0 O - e f (£> q) 


(62a) 

(62b) 
The \ f s 

(63a) 

(63b) 

(64) 

(65a) 

(65b) 

( 66 ) 
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where: 


Va + b 


q 


2 


2b 

a + b 



( 67 ) 


and f(;J, q) is the incomplete elliptic integral of the first kind. If b > a > 0, Equation 288.50 of Refer- 
ence 2 gives: 


2 V 0 1 gf(0, q) 


( 68 ) 


where: 


Define: 



a f b 

^ ~2tT 


^ sin_1 /a+ b 


k 90° - T 


(69) 


then: 




J .rr/2 

o <«*'-»> 


(IT 

cos T) 1 


(70) 


By the last of Equations (64), a' >b. For a' >b> 0, Equation 291.00 of Reference 2 gives: 


2 ') 0 ^ gf(A, q) 


(71) 
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where: 


g s 


ys^Tb 


2b 

a ' + b 


A = 


1 / 

Vq 2 ( 


a ' ~ b) 


The incomplete elliptic integral of the first kind is extensively tabulated. 
Finally, the form in Equation (27e) was: 

A' + B' cos 2L + C' cos 4L 


( 72 ) 


In our particular case: 


A' 


B' = E 


1 - E 2 {\~V +| 

-\v) 




we may write: 

(A ' + B ' cos 2L + C ' cos 4L) ~ 1/2 - (l + E-^ 2 Ecos 2 l) 1 2 (1-E + / £ 2 Ecos 2 l) 1 


\ 0 + 2k 2 cos 2L + 2 \ 4 cos 4L + 


write: 


(1+E-^ 2 Ecos 2 l)” 1/2 - 77 0 - 2t] 2 cos 2L + 2r^ cos 4L - * * * 


(l - E + 't 2 E cos 2 L) 


1/2 


6 o + 26 l 2 cos2L + 26 ) 4 cos 4L + 


(73) 


(74) 


(75) 


The tj’s and 0’s are related to the Vs as in Equation (46). The second factor of Equation (74) may 
be written as: 


(l -E + 't 2 Ecos 2 L) = (l - E + t 2 e) - j_g + E ^ sir > 2 ^ 


(76) 
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which is in the form of Equation (28). The first factor of Equation (74) is: 


(l + E - 'C 2 E cos 2 l) 


(1 + E) 


( 2 \ 

II - YJx cos 2 L J 


Define: 


L * 90° - T 


Then: 


(l +E -l 1 Ecos 2 L) 


(1+E) 



-t 2 E 

1 + E 


and this is in the proper form. 
From Equations (75): 




e 


0 


2 

77 


2 

77 


J * 77-/2 
0 



dL 

(l + E -'C 2 Ecos 2 L) 1/2 


dL 

(l - E + -£ 2 E cos 2 L ) 1/2 


( 77 ) 


(78) 


(79) 


Equation (78) is transformed to the form of Equation (28) by the substitution, L = 90 - T. Equation 
(79) is transformed by simply writing cos 2 l = l-sin 2 L. 

We have seen how the partial anomalies are defined and how the various radicals can be ex- 
panded in a form convenient for integration. The partial anomalies will permit us to obtain highly 
convergent series by means of which the perturbations of a body may be expressed. The require- 
ment that the series be highly convergent will determine how we must select the points of separa- 
tion in dividing the ellipse. 

The series for computing perturbations, when expressed in the conventional anomalies, may 
become very slowly convergent or divergent if the radius vector of the perturbed body is at one time 
much larger and at another time much smaller than the radius vector of the perturbing body. Con- 
sequently, the points of separation must be selected such that in each segment of the perturbed 
ellipse, the radius vector of the perturbed body is always larger or always smaller than that of the 
perturbing body. If the two bodies approach each other closely, the small segment of the perturbed 
orbit near this minimum of distance should be represented by the intermediate anomaly. If the 
segment represented by the intermediate anomlay is small, the value of the modulus, l, will be 
very small and this will lead to highly convergent series. 
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m. DETERMINATION OF CONSTANTS OF INTEGRATION 


We now consider the types of integrations we will be performing and the determination of the 
constants of integration. As we have seen in Section n, the radicals involved in the various partial 
anomalies may be expressed in the general form: 


sin k + J3 2 cos 2k + /3^ sin 3k + (3^ cos 4k + * * * (80) 

In certain cases the coefficients of the sine terms are zero. The expressions for the polar coor- 
dinates of a body in an elliptic orbit have been given in Section I in terms of the partial anomalies 
and these coordinates may also be expressed in the form of Equation (80). Since the perturbation 
function and its derivatives involve these coordinates, they may also be expressed in this form. 

Assume that the perturbed body is a comet and the perturbing body a planet. Let g' be the 
mean anomaly of the planet. If y is any elliptic element, we may write: 


dy 


i= o 


+ a. 


sin k + a 2 cos 2k + 


\ sin 
' cos 


ig 


(81) 


The symbol ig' implies that a series in sin ig ' applies when y is a particular elliptic element 
representable by a sine series, nor a>, for example, and a series in cos ig' applies when y is a par- 
ticular elliptic element representable by a cosine series, a, e, i, for example. 

The results of Sections I and II show that n dt may be written as: 


n dt - (\j cos k + k 2 sin 2k + \ 3 cos 3k + • ■ • ) dk (82) 

where n is the mean motion of the comet defined by Kepler's third law. Eliminating dt between 
Equations (81) and (82): 


dy 



cos k + * 2 sin 2k + • • • ) 


cos . 
sin 


(83) 


If n’ is the mean motion of the planet and c* is its mean anomaly at t =0, we have: 

g' = n' t + c' (84) 


Define: 


V 


n 


(85) 
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Integrating Equation (82) and using Equation (85), we have: 


n't - vk Q + sink - 2 2 cos 2k + 3 zA 3 sin 3 k ^ * * * (86) 

where k Q is the constant of integration. We now eliminate g' from Equation (83) by means of 
Equation (86) and obtain: 


d Y r dk { K i 

4*0 


cosk + k 2 sin 2k + k 3 cos 3k + 


5in /. , 

:os i lc 


+ i i\ Q + iiAj sink - y izA 2 cos 2k + 


•) (87) 


By known trigonometric identities: 


sin (ic ' + izA Q + * * * ) " sin ic' (f9 0 + 9 x sin k + 0 2 cos 2k + • * * ) + cos ic 1 ( Oq +6^ sin k + 0 2 cos 2k + * * * ) 
cos (ic ' + izA 0 + * * ’ ) ~ cos ic ' (# 0 + 0 x sin k + 9 2 cos 2k + • • • ) - sin ic 1 (d Q ' + 6 ^ sin k + 0 2 cos 2k + * * * ) 

Equation (87) becomes: 


dy - dk n cos k + ct > 2 sin 2k + <<? 3 cos 3k + • ■ ■ ) 


( 88 ) 


Since c' is independent of the time and hence of the partial anomaly, we have: 


Y 


const 



sink 


1 

2 u 2 cos + 


1 

3 > 3 s in 3k 



sin 

cos 


(89) 


Let us select a particular perihelion passage of the comet as the epoch of time. Then, g' = c' 
at t = 0, and c 0 ' is the mean anomaly of the planet at the instant of this particular perihelion pas- 
sage of the comet. Let T be the time of the closest perihelion passage of the comet, past or future, 
from the time t. Then, c' will be the mean anomaly of the planet at time T. Thus, c' is fixed 
during each revolution of the comet but changes from revolution to revolution. 

In order to proceed, we must select the points of separation. To illustrate a fairly complex 
example, we select six points of separation. 

These points are: 


r - a( 1 - e) 


r 
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1 


r - a( 1 + e ) 
r = r" 


This is illustrated in Figure I. We have selected f' and f / less than 180° and f" and f 
than 180°. The following partial anomalies are required. 


Segment of Ellipse 
Perihelion to r ' 

r 1 to r/ 
r j' to Aphelion 
Aphelion to r" 
r" to rj" 

r/' to Perihelion 

We write Equation (89) for each segment. 


Partial Anomaly 

K i 

L i 

k 2 

k 3 

k 4 


Y - V®) + 2>< 2) cos2K 1+ ,^cos4K 1+ ...)^ic' 

Y = Xj (0) + cos 2L, + cos 4L, + • • •) ^ ic' 

Y = « 2 (0) + ^(' < 2 (2) cos 2K 2 +K 2 (4) cos4K 2 + ‘ ") sin ic ' 

Y = * 3 (0) + 2Z( K3<2> C ° S 2K} + * 3<4> Cos4K 3 + •") sin ic ' 

Y = ^4 (0) + 2^( X 4 (2> cos 2L 4 +X 4 <4) cos 4L 4 + • • •) sin ic ' 

Y = K 4 <°> + 22 (« 4 <2) cos 2K 4 + cos 4K 4 + ■ • •) ^ ic* 


greater 


(90a) 
(90b) 
(90c) 
(90d) 
(90e) 
(9 Of) 
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K i°\ k i 0) 9 k 2 (0) > k 3 <0) > k 4 0) > * 4 (0) are constants of integration. We can see from the results of 

Section I that for this division of the ellipse no sine terms are present. 

Define: 


+ S 1 

- 


+ 

+ k/ 6 ) 4 ••• 

(91a) 

+ s 2 

= 


+ K 2 («) 

+ k< 6 ~> + 

(91b) 

+ s 3 


«3 <2) 

+ kW 

+ k 3 <«) 4 

(91c) 

+ s 4 

= 

S (2) 

4 K < 4) 

4 

+ u (6) + ••• 

(91d) 

+ Ui 

- 


+ 

4 \W 4 ... 

(91e) 

+ u 4 


\ <a> 

4 

+ \ (4 > 

4 \ <•> 4 * * * 

4 

(91f) 

- S 1 


K<~V 

- X W 

+ K( <6) .... 

(91g) 

- s 2 



- x 2 «> 

4 *,<•> - 

(91h) 

- S 3 

= 


- 

+ „«> - ... 

(91i) 

- s 4 

- 


~ X 

*4 

+ k 4 (6) 

(91j) 

- U 1 

= 


- ^< 4 > 

4 

(91k) 

- 

= 

v 2) 

- JV < 4) 

4 

4 X (•> - • • • 
4 

(911) 


Let us select a particular value of c', say: 


c - c r 


At the epoch corresponding to c 0 ', K 2 = 0° since we are at perihelion. Equation (91a) permits us 
to write Equation (90a) as: 


Y ~ K 


( 0 ) + V 1 M cos 

1 L — J sin 


W - (Y) 
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In this manner, (y) is defined as the value of y at t = T 0 . Thus: 


( 0 ) = 


(Y) 


\ ’ COS 

/ . M i sin ic o 


( 92 ) 


At r - r' , Kj - 90° and Equation (90a) by means of Equation (91g) becomes: 


Y r ' 


«/°) 


\ ’ COS 


But at r = r', Lj =0° and Equations (90b) and (91e) give: 


V = \ (0) + 2I P isin ic o' 


Since these expressions must be identical, we have: 




Z>.: 


( 0 ) 


-& 


Combining this with Equation (92): 




(Y) 



■ p .) 


COS . 

sin ic o 


(93) 


We proceed in this manner at the other points of separation and obtain the following equations 
for the constants of integration. 


*i (0) = (Y) -£ M lsin ic o 

\ (0) = (Y) 'Z]K +N i + Pi)sin^o' 

k 2 W " (Y) +N 1 +P 1 +< ?i" N 2)sin ic o' 

= (Y) -2Z( M ^ +N 4 +P 4 +Q 4- N 3)sin ic o' 


(94a) 

(94b) 

(94c) 

(94d) 
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(94e) 


\ (0) = (Y) - 2I( M 4 +N 4 +P 4)sin ic 0 

K 4 < 0) = (Y)-^M 4 ^ic 0 ' (94f) 

During the x th revolution after T 0 , we have: 

c . c; 

<Y) X = Y at t = T X 

Then in the same manner as Equations (94) were obtained: 


k/ 0 ) 

S * cos 

= (Y), ~ / M 1 sin 1C , 

(95a) 

K (0) 

= (Y), "2j M l +N . +P l) stnic,' 

(95b) 

K <°> 
K 2 

= (Y) x 'Z]( M l +N l +P * +Q l ' N 2 )sin ic x' 

(95c) 


= (Y) x -Z]( M 4 +N 4 +P 4 + Q4- N 3)stn iC X ‘ 

(95d) 

k <°> 

4 

= (Y), -2J M 4 +N 4 +P 4)sin ic / 

(95e) 

K (0) 

K 4 

= (Y) x -2I( M 4)sin ic ,' 

(95f) 


We must now relate ( y ) , (y),, • • • (y) x . Consider the first aphelion passage of the comet 
after t = T 0 . Using Equation (94c) for /< 2 (0) in Equation (90c) gives: 

Y = (Y) +N X +P, +Qj -N, -M 2 ) sin ic 0 

= (Y) - 2^( Ri +Tj -R 2 ) ^ ic 0 ' 
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where we have recognized that at aphelion, k 2 = 180°. But, at aphelion k 3 = 180° and we may con- 
sider this to be the beginning of the revolution with c' = c,' . Equation (95d) gives « 3 (0) with x = 1 
and using this value in Equation (90d), we have: 

Y = (Y)! - 2^](R 4 +T 4 -R 3 )^ic 1 ' 

Since these two expressions for Y must be identical, we have: 


(Y), = (Y) - 2^(R, + T, -R a )“*ic 0 ' 

- 2 ZX +T 4- R 3)s?n ic l* 

In general, for any two consecutive values, (y) x _, and (y) x , we have: 

(Y) x = (Y)^ - 2 2](R 1 + T 1 -R 2 )^ic x '. 1 

‘ 2 + T “ ~ ^in ic x (96) 


We use Equations (96) to relate (y) to (y) and find: 


(Y) x = (Y) 


- 2^(8! + T 1 -R 2 ) gi 


COS . cos . 1 

lc o + sin 1C l + **’ + sin ic x-ij 


\ ^ \ r cos . , cos . , cos *1 

2 / J R 4 + ^4 R 3 ) L s in 1C l + sin ic 2 + * * ’ + sin ic x J 


(97) 


Let A be the mean motion of the planet during a mean revolution of the comet. Then: 


c 0 # + A 


c/ + A 


c o + 2A 


c ' = c' + xA 

x 0 


(98) 


i 
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Two trigonometric identities are useful. 


1 + cos iA + cos 2 iA + • * * + cos (x - 1 ) iA - y + 


1 s in (x - 2 ) iA 


2 sin j iA 


;in iA + sin 2 iA + * * * + sin (x - 1 ) iA - 2 cot y iA 


(x- j)iZ 


2 s i n 2 i A 


Multiplying the first identity by cos ic 0 ' and the second by sin ic 0 ' and subtracting, we have: 


cos ic Q ' + cos iCj' + ** * + cos ic^ - y cos ic Q ' - 2 cot A sin ic Q ' + 2 cot ^ A sin ic x ' - *2 cos 1C x 
Multiplying the first identity by sin ic 0 ' and the second by cos ic 0 # and adding: 

. , . . , . . , 1 . . , 1 i . , 1 i . ( 1 . , 

sin ic Q + sin iCj + *■* + sin ic H - 2 sin ic Q + ^ cot ~2 ^ cos lc o ~ 2 cot 2 A cos ic x - | sin ic x 

These expressions are to be substituted in Equation (97). We must separate the sine and cosine 
terms. The quantities with subscript c are to be used when y is given by a cosine series and the 
quantities with subscript s are to be used when y is given by a sine series. Equation (97) becomes: 


(*) x = (Y) -Z]( R «c +T .c- R 2 c +R 4c +T 4c- R 3c) COsiC 0' 

+ l>ic +T lc _R 2 c _R 4c _ T 4 c +R 3c) COt 1 Asin iC 0 


-/ (R 1 + T. -R, -iR +T -R ) sin ic n # 

/ 1 \ls Is 2s 4s 4s 3s/ 0 

" / (R, +T. “R- -R -T +R ) cot 7 T A cos ic ' 

/ vis Is 2s 4s 4s 3s/ 1 0 


+ TK +T lc " R 2c +R 4c +T 4c _R 3c) COS iC x' 

" +T > " ‘ T 4c +R 3,) COt 1 Asin iC x 

+ +T <s - R 2 x - R 4 s- T 4s +R 3 s )^n i ‘-'x‘ 

4 2Z^ R >» +T ‘ s " R 2* " R ^ ~ T 4 s 4R 3 S ) COt 1 AcOS > C x 
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As c 0 ' is a constant depending on the initial conditions, its effect is included in (y) which must 
be found from astronomical observations. We omit all terms dependent on c 0 ' . The subscript x 
on c' is dropped and c' is the mean anomaly of the planet at the perihelion passage of the comet 
closest to the time at which we are determining the perturbed orbital elements. 

Taking account of the above remarks and using Equation (99) for (y) x , Equations (95) become: 

^ (0) S = (Y) + 2^ S l C +T lc- R 2c +R 4c +T 4c- R 3c + X s ) C °S 

+ S ls +T ls " R 2s +R 4s +T 4s _R 3s _X c) sin ic ' (100a) 

V 0> = (Y) + 2j- R lc - U l. - R 2 c + R 4 c + T 4c - R 3 c + X s ) COS ic< 

+ R ls ~ U 1 S “ R 2 s +R 4s +T 4s " R 3s ~ X c) Sin iC ' (100b) 

"2 (0) = (Y) + Z](" R lc- T lc- S 2c +R 4c +T 4c- R 3c +X s) COsic ' 

+ " Tls “ S 2 s +R 4s +T 4 s ' R 3s * X c) Sin ic ' (100c) 

x A (0) = (Y) + ^( R lc + T lc -R 2c -S 4c+ T 4c -R 3c+ X s )co S ic' 

+ +T 1 s +T U ~ R 3 s ~~ X c ) sinic' (100d) 

V 0) = <Y) + 2^( R lc +T lc- R 2c- R 4c- U 4c- R 3c +X s) COSiC ' 

+ +Tls “ R2s " r «» ~ U 4 s ~ R3s " X ^ sinic ' (100e) 


4 (0) = (Y) + ^(R lc + T lc -R 2c -R 4c -T 4c -S 3c+ X s )cosic' 


Ls _R 2s " R 4s _ T 4 s ~ S 3s " X c ) sinic ' 
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where: 


X c " ( R lc +T lc - R 2c - R 4c - T 4c +R 3c) COt T A 
X s S ( R ls +T ls ~ R 2 S ~ R 4s ~ T 4 s +R 3s) COt 1 A 

In general then: 

Y " (Y) + i (0> + *i <2) cos 2K + * * * ) sin ic' 


If for some value of i, say i = i f , we have: 


i f A -- m • 360° 


where m is an integer, cot (i/2) A becomes infinite and the above expressions are useless, 
ing to Equation (98) together with Equation (103): 


= c o + r 


360° 


This may be written as: 


if C X = if c 0 ' + mx • 360° 


Since mx = 0, 1, 2 • . we have: 


Equation (97) becomes: 


(Y) x 


(y) - 2x(R lc + T lc -R 2c -R 4c -T 4 , 


- 2x ( R . s +T Is - R 2s* R 4,- T 4 


+ R 3l .) cos i f c' 
+ r 3 J sin i f o' 


(101a) 

(101b) 

( 102 ) 

(103) 

Accord- 


(104) 

(105) 
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Equations (95) become: 



K< °) 

- ( Y ) “ ( R lc +S lc) cos i f c ' ■ xz c cos 1 f c ' 




~ ( R ls +S ls) sin if c ' - xz s sin i f c ' 

(106a) 


\j(0) 

= (Y) ' (2R lc 

+ T. + U, ^ cos i , c ' - xz cos i , c 1 

1 c 1 c / f c f 




' (2 R ls 

+ T. + U. ^ sin i , c ' - xz sin i t c' 

lsls/ f s f 

(106b) 


k 2 <°> 

= (Y) - (2 R lc 

+ 2T lc -R 2c +S 2c ) cos i f c' - xz c cos i f c' 




- (2 R ls 

+ 2T. -R 9 + S 0 ) sini c' - xz sin i f c' 

Is 2s 2s/ f s f 

(106c) 


«3 <0) 

= (Y) - (2R 4c 

+ 2T d - R, + S, ) cos i . c ' - xz cos i , c ' 

4c 3c 3c/ f c f 




- K. 

+ 2T d - R, +S, ) sin i f c ' - xz sini. c' 

4s 3s3s/ f s f 

(106d) 


\<°> 

= (Y) - (2R 4r 

^ T„ + U „ ) cos i f c' - xz cos i , c' 

4 c 4 c ) f c f 




" (2R 4s 

+ T. + U , 1 sin i , c ' - xz sin i ( c' 
4 s 4 s ) f s f 

(106e) 


* 4 <°> 

= (Y) - ( r 4c + S 4c ) cos i f c ' " xz c cos i f c' 


where: 


- ( r 4s +S 4s )sini f c' - xz s sini fC ' 

(106f) 



z c s 2 ( 

Rjc +T lc - R 2 c " R 4c "" T 4 c +R 3c) 

(107a) 



z s - 2( 

»!, +T ls - R 2 S - R 4 s - T 4s + R 3s) 

(107b) 


This special case will always arise for i f = 0. 

In order to calculate the perturbations of the elliptic elements involving the mean motion, we 
must integrate again. If y represents the first integral obtained above, the integral, J y dt is re- 
quired. Let us call this integral z. 


37 


The general form of y is : 


Y 


(Y) 


2>“’ + h* ^ sink + h^ 2 ^ cos 2k + * * * ) 


(108) 


The sine terms may be absent in some cases as they are for our division of the ellipse into six 
segments. The form of dt is given by Equation (82) and hence y dt is of the form: 


The integral is: 


dk 2j^ cos k + 4> 2 s in 2k + * * * ) 


const + 


w 


sink - 2 tf> 2 


cos 2k + 


> cos 
' sin 


(109) 


This is exactly the form of the first integral, Equation (89). The constants of integration can be 
found in the same manner as for the first integral. The special case, i f A 2™, requires separate 
consideration with regard to the terms multiplied by x. The terms independent of x are dealt with 
as in the first integration. The terms in x are of the form const x x, and dt is of the form of Equa- 
tion (82). Since x is independent of the partial anomaly, we will have the form: 


z 



const + x[/dj cos 2K + • • •] 


where we are considering only the terms involving x and, since in our case sine terms do not ap- 
pear, only cosine terms are included. We will have an integral of this form for each segment of 
the ellipse. Thus: 


z _ j + x (j cos 2Kj + j t (4) cos 4Kj + • • • ) 

z - ■tf 0 ^ + x(-t,< 2 > cos 2Lj + cos 4L, + • • •) 

z - j 2 (0) + x (j 2 (2) cos 2K 2 j 2 (4) cos 4K 2 + * ‘ •) 

z - j 3 (0 > + x(j 3 <2> cos 2K 3 + j 3 (4> cos 4K 3 + • • •) 

z 4 x(t’ 4 <2) cos 2L 4 + f 4 <4> cos 4L 4 < • • •) 

z j 4 (0) 4 X (i 4 < 2 > cos 2K 4 + j 4 (4) cos 4K 4 ■* • • • ) 


(110a) 

(110b) 

(110c) 

(llOd) 

(llOe) 

(llOf) 
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As in Equations (91), define; 



= r; 

+ s,‘ 

= ii (2) 

+ j , (4) 

+ • • • 

* 


m 2 ' 

= r 2 ' 

+ S 2 ' 

= h (2) 

+ j 2 <4) 

+ • • • 

, etc. , 


N,' 

= .r; 

- s ,' 

= > r 

- i< 4 > 

+ • • • 

■ 


n 2 ' 

= r 2 ' 

- S 2 ! 

J 2 

- 


, etc. 

(Ill) 


Let (z) x be the value of z at the x th passage of the perturbed body through its perihelion after 
the epoch. A calculation analogous to that which led to Equations (94) gives: 


= (Z) x - xM; (112a) 

V 0) = (*), - x(m,' + n; +Pl ') (112b) 

i 2 <0) = (z) x - x(m; +n; +p/ + Q/ -n 2 ') (H2c) 

if = ( 2 ) x - x(M 4 * +N 4 ' +p; +Q;-N 3 ') (112d) 

V 0) = (z) x - x(M 4 ‘ + n; + P 4 ' ) (112e) 

i 4 (0) = (z) x -xM 4 ' (112f) 

Also: 


(Z) x = ( Z )x-1 " 2(x- l)(Rj' +T,‘ -R 2 ') + 2x(r; +T 4 * -R/) 

which gives: 

(z) x - (z) - 2 [(x ~ 1) + (x - 2) + * * * +2 + l](R/ +T; *R 2 # ) + 2 [x+(x-1)+---+2+i](r; + T/ - R 3 ')(113) 

But: 

1 , l 

l + 2 + -- '+ (x-l) - ~2 x 2 - "2 x 
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and: 


1 + 2 + • ' • + x - 2 x2+ 2 


Equation (113) becomes: 


Define: 


Equations (112) become: 


- x 2 (r; 

+ t j* - r 2 -r 4 '-t; +r 3 ') + x(r; +t,'-r 2 ' +r; +t; -r 3 ') 

(114) 


w = R/ + T/ -r 2 * -r; -t 4 ' + r 3 ' 

(115) 

j/o = 

(z) + x(- S ' + T,' - R 2 +R 4 ' + t; -R 3 ') - x 2 w 

(116a) 

iw = 

(z) + x(-r; -u,* -r 2 ' +r 4 ' + t 4 ‘ -r 3 ‘) - x 2 w 

(116b) 

j 2 (0 > = 

( Z ) 4 x (-r; -Tj'-Sj' +r; + t; -r 3 ') - x 2 w 

(116c) 

j 3 <0> = 

(z) 4 x (r 1 ' 4 t; -r 2 '-r; -t;-s 3 ') - x 2 w 

(116d) 

<<°> - 

4 

(z) 4 x (r 1 * +t; -r 2 '-r; -o 4 '-r 3 ') - x 2 w 

(116e) 

j (°) = 

J 4 

(z) 4 x(R l ' 4 T ; -r 2 - -s; +t; -r 3 ') - x 2 w 

(116f) 


The determination of the constants of integration has been carried out for a division of the 
ellipse into six sections using the partial anomalies K,, K 2 , K 3 , K 4 , L,, and L 4 . It is of interest to 
determine what changes occur for different divisions of the ellipse and for the use of different 
anomalies. 

If we use L 2 and L 3 rather than L, and l 4> we must replace: 

P,. Q,. P 4 - T r T 4* U l’ U 4 


by: 


-q 2 , -p 2 . -q 3 . -p 3 . -t 2 . -t 3 . + u 2 , ^3 


respectively. 
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If we used Lj and L 2 instead of L x andL 4 , we would make the above changes only for those 
constants relating to L 4 . If we used more intermediate anomalies to divide the ellipse into more 
than six parts or more than two parts using these anomalies, we simply increase the number of 
p T s, Q's, etc. 

If we omit certain points of separation, we simply omit the corresponding constants unless we 
omit the point of separation at perihelion or aphelion. If we use the same partial anomaly on both 
sides of perihelion or aphelion, i.e., if we have no point of separation at perihelion or aphelion, the 
equation representing the element y will be of the form of Equation (108), that is, sine terms will 
be present. This will change the form of the equations. 

Suppose, in our example, we omit the point of separation at perihelion. In all of the equations 
we will haveMj = M 4 . In addition, Equations (90a) and (90f) will combine and take the form: 

y = ^ sink+/< (2 ) cos 2k + **’) s j n ic # (117) 


Define: 


E = G + H 


*(D + *<2) - *(3> _ k (4) + ... 


(118a) 


F = G - H 


K d) - „<2) _ *<3> + k (4) + «(5> 


(118b) 


Let us omit Equation (90f) and all quantities dependent on it. Equation (90a) in the form of 
Equation (117) describes this segment. We see that in all of the equations for the constants of 
integration, we drop all terms with subscript 4 and change: 


to: 


M r N r R r Sj 

*<°\ -E, -F, -G, -H 


Now, omit the point of separation at aphelion. Equations (90c) and (90d) will combine and take 
the form: 


Define: 


\ i cos 

Y = kW + / >,<»> sink, + « t (2) cos 2k, +••■) sin ic* ( 119 ) 

E, s G, + H, - kW + - k/ 3) - + ... (120a) 

F, 3 G, ~ H, 3 «, (l) - «j (2) - /</« + k/ 4 ) + k/ 5 > - ••• ( 120 b) 
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Let us omit Equation (90d), and consequently, drop all terms with subscript 3. Then we must 
change: « 2 (0> , m 2 , n 2 , R 2 , s 2 to: 


S (0) . "El- -F,. -0.. -Hi 

In this case, we must take x to be equal to the number of complete revolutions of the comet since 
the epoch if the time we are considering corresponds to the anomalies k 1? L t , or , and must take 
x to be equal to this number plus one if the time we are considering pertains to L 4 or k 4 . 

We could also omit Equation (90c) and write Equation (90d) as: 

Y = K n 0)t rH" Sink l +K l ( l 2> cos2k l + "') sin ic ' • (121) 


Define: 


E U r G 11 + H U K 11 1> + K H 2) - *n 3> - *1 ( 1 4) + ••• (122a) 

F„ = On " H n " k/, 0 - «{*> ~ « t \ 3) * «, { i 4) + ••• (122b) 

consequently, we drop all terms with subscript 2 and change: * 3 (0) , M 3 , n 3 , R 3 , S 3 to: 

K (°) F F G - H 
11 ’ r ll ’ ^11 ’ u ll 1 n 

In this case, x has its usual definition when using K 1 and L t , but must be increased by unity when 
using k j, L 4 , andK 4 . 

In the method of partial anomalies, the integrations have not led to the divisors (in - i' n') and 
(in - i* n' ) 2 which appear in other methods. Instead, we obtain the multipliers cot (i/2) A and 
cot 2 (i/2) A. In the case of near commensurability, these multipliers become large and lead to 
large perturbations of long period. When in - i ' n' is small, it may be that both i and i' have large 
values and, consequently, by ordinary methods we must develope many terms of the series in both 
i and i' . In the method of mean anomalies, this is necessary for only one of the indices because 
of the rapid convergence of the series involving the partial anomalies. 
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IV. THE PERTURBATION EQUATIONS 


In the following, the variables y, H, V, p, and q, first introduced by Hansen (Reference 3), will 
be used to designate the elliptic elements. 

The perturbation equations are: 


5k 

dt 


2 an / a 1 \ 

VT=? \ r + l-e*) 


c dQ 

cosf + 


e an 

1 - e 2 


+ y sin f 



(123a) 


dW 

dt 




sin f 


an a 

If - 7 cos f ' 


an 

dr 


(123b) 


dH 

dt 


3 an dH 


(123c) 


dp 

dt 



en 


cos i sin I sin f sin ( f ' + v - k) ^ 


(123d) 


dq 

dt 




cos i sin I cos f sin ( f ' + v - k) 


cM2 

<9H 


(123e) 


In addition: 


dSz 

dt 


Y ~ cos f + ” sin f + H 


(124a) 


dw 

dt 


2 n - 7 


sin f - m 


2 /r 


(cos f + e) 


(124b) 


Ss D Sq sin f - Sp cos f 


(124c) 


In these equations, nSz is the perturbation in the mean anomaly or mean longitude, w is the per- 
turbation in the natural logarithm of the radius vector, and 5s is the perturbation in the sine of 
the latitude. The elements y, and H are obtained by integrating Equations (123a), (123b) and 
(123c), and Sp and Sq are the integrals of Equations (123d) and (123e). Also, I designates the in- 
clination of the orbit of the comet with respect to the orbit of the planet, i is the inclination of the 
orbit of the comet with respect to the ecliptic, and 5s is measured with respect to the ecliptic. 
The quantity ft is the perturbation function defined by: 


ft = 



(125) 
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m' and r' being the mass and radius vector, respectively, of the perturbing body, H the cosine of 
the angle between the radius vector r and r ' , and A is the distance between the planet and the comet. 
In addition, v + k is the angle measured in the plane of the orbit of the comet, and measured from 
the ascending node of the orbit of the comet on the orbit of the planet to the perihelion of the comet, 
while v -k is the angle measured in the plane of the orbit of the planet, and measured from the as- 
cending node of the orbit of the comet on the orbit of the planet to the perihelion of the planet. 
Consequently: 


H ^ cos(f + y + k)cos(f'+^-k) + cos I s in ( f + ^+ k)sin(f' + u ~ k) 


( 126 ) 


Also: 


A 2 = r 2 + r ' 2 - 2rr' H 


( 127 ) 


We also have: 


0 - 3r cos f dy + 3r sin f dM* + 4a dH ( 128 ) 

w = - | y | cos f - j V “j- sin f - Js ( 129 ) 


write: 


n - Qj + n 2 


Then: 




(Jfij 


r r ' 

dH 


Tf 

m ' 

A 3 

d f 


dCl l 


rr ' 


r 2 

r TT 

- m ' 

A 3 

H - m' 

A 3 

an x 


r r ' 




- m' 

A 3 




( 130 ) 


(131) 
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By Equation (126): 


dH 

df 


sin (f + v +k) cos (f ' + v -k) + cos I cos ( f + v + k) sin ( f ' + v - k) 


(132) 


Equations (131) become: 


<9Qj 

r r ' 


- ~m' , 

A 3 

an, 

r r ' 

d r 

- m ' 

A 3 

dQ, 

rr ' 

~dW 

m ' “ 

A 3 


|sin(f + i' + k)cos(f , +z^-k) - cos I cos (f +j/+k)sin(f' + V - k)J 

r r ' r r 

A 3 L COS ( f + 1 ' + k ) cos ( f ' +l/_k ) + cos I sin (f + v + k) sin ( f + v ~k)J -m' “ 


2 

A 3 


(133) 


and: 


Define: 

cos ( V - k) = y sin G 
cos I sin (i/- k) = y cos G 

v + k + G - r 

- sin (v - k) s y' sin G' 

cos I cos (y - k) = y'cosG' 

v 1 k f G' = F' 

Let u ' be the eccentric anomaly of the planet. The elliptic equations give: 

r ' cos f ' ~ a ' cos u ' - e ' a ' 
r' sinf' = a' Vl ” e ' 2 sinu' 


(134a) 


(134b) 


(135) 


Equations (134), and (135) allow Equations (133) to be written as: 


~JT 



y cos ( f + T) cos u / - e ' y cos ( f + F) 
+ y' /l - e ' 2 cos ( f + F' ) sin u ' 


(136a) 
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r 


a r 


7 sin ( f + P) cos u' - e ' 7 sin ( f + T) 


(136b) 


<?Qj 

~dT 


A 3 


+ 7' VT-e ' 2 sin ( f + T' ) sin u ' 


r 

a ' 


_ ' rr ' 

dH " m ' A 5 


Also: 


r'sin(f'+v-k) - -a' 7' sin G' cos u' + e' a' 7' sin G' + a' 7 Tl- e' 2 sin G sin u ' 


Equations (123) may then be written: 


where: 


dv £ 

dk = m ' 


A 3 


|^2 P cos u' + 2Qsinu' + Rj 


d*y 

dk “ m 


' ~ ^ 2 Pj cos u ' + 2 Q t sin u ' + R J 


dk " m ' ^3 [ 2 P 2 cosu' + 2 Q 2 sinu' +R 2 ] 


" m ' ^3 [ 2P 3 cosu ' + 2 Q 3 sinu' +R 3 ] cos i 


dq _ r 3 r -j 

dk - m' — |^ 2 P 4 cos u' ^ 2Q 4 sinu' +R 4 j cos 1 


a ' a 2 7 ndt 

r Vl ~ r 2 dk 


cosr + ^7 >s f ’- p cos ( f + r> 

a(l-e 2 ) 


! 7 Vl - e ' 2 ndt 


r 2 dk 


r(cos f + e) 

cos F' + Tv — cos ( f + f ) 

a (l - e 2 ) 


2a ' a 2 7 ndt 


Vl ~ e 2 r 2 dk 


( cc 


. r(cosf + e) \ r sin f 

e,COSr + “ a (i“e'i)“ cos ( f +l >) + - 77 “ 


/ 2 

a a 7 


ndt 


r Vl - c 2 r 2 dk 


r sin f 

sin P - — ” " cos ( f f T) 


(1- 2 ) 


(136c) 

(137) 

(138a) 

(138b) 

(138c) 

(138d) 

(138e) 

(139a) 

(139b) 

(139c) 

(139d) 
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Qi 


P 3 ° 


a' a 2 y' Vl-e' 2 

ndt 

r Vl -e 2 

r 2 dk 

2a ' a 2 y ndt 

e ^sin r 

r i\ - e 2 r 2 dk 

3a' a y ndt 

‘ iVT^ ’ r 2 dk C ° S(f + 

3a ' a y* VT- e ' 2 

ndt 

2 i\ - e 2 

r 2 dk 

3a' aye' ndt 

/l “ e 2 r 2 dk 

cos ( f + r 

a ' a y* sin G ' sin I 

ndt 

2 i\ - e 2 

r 2 dk 

a' ay V 1 - e' 2 sinG sin I 

2 Vl - e 2 

e * a ' ay' sin G ' 

sin I 

Vl “ e 2 

r 

a ' ay' sin G ' sin I 

ndt 


r sin f 

Sin r - — 2 \ cos ( f + r # ) 


1 (l - e 2 ) 


r sin f 


a 


(1 - e 2 ) 


cos ( f + r) ) + 


sin f 


r 2 dk 


sin f 


sin f 


r 2 dk 


cos f 


r cos f 

a ' y 


^ a' ay j \ - e ' 2 sin G sin I ndt 

V 4 „ r o * „ .. COS f 


2 /l - e 2 


r 2 dk 


e ' a ' ay' sin G' sin I ndt 




r 2 dk 


cos f 


Expressions for ndt in terms of the various partial anomalies were given in section I. 
We must now take account of fi 2 given by Equation (130). Recall: 


fi, = -m' — y H 

r 


(139e) 

(139f) 

(139g) 
(139h) 
(139i) 
(139j ) 
(139k) 
(1391) 
(139m) 
(139n) 
(139o) 


47 


Then: 


_ m' r dH _ m ' rr ' <9H 
d f ~ _ r ' 2 d f ” ~ r > 3 <9 f 


(140a) 


dCt 2 

r ~J7 



m' rr ' 


r 


' 3 


H 


(140b) 


dQ 2 



(140c) 


If, in Equations (131), we change A 3 to ~r' 3 and omit the term-m' r 2 /A 3 ? these equations are iden- 
tical to Equations (140). The results for the second part of the perturbing function are the same 
as for the first part if in Equations (139c) and (139f) we omit the term r sin f 'a' y and r cos f/a' y , 
respectively, and if, in all equations, we replace A 3 by - r' 3 . Consequently, Equations (138) and 
(139) show that the expressions for dy/dk, dy/dk, etc., arising from the second part of the perturb- 
ing function, contain terms of the form: 


CO S U ' - c ' 


r 


/ 3 


, sin u 

and — — 

r J 


These expressions may be developed in a Fourier series in cos iu ' or sin in'. The coefficients of 
these series will be given by: 




1 [ ( co s u ' - s i n <£' ) cos i u ' du ' 

a' 3 rr J Q ( 1 - sin 0' cos u' ) 3 


1 I sin u ' sin iu ' du ' 
a>3jT J 0 (1 ” sin 0' cos u' ) 3 


(141a) 


(141b) 


where: 


e ' - sin 0' 


(142) 


and the elliptic equation relating r ' to a', e' , and u' has been taken into account. Integration by 
parts gives: 

r v 

1 I cos iu ' du ' 

n J 0 1 ~ sin 0' cos u 1 


1 , 
tan 1 ry 

— (143a) 

COS(/)' 
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Differentiating this with respect to 0 ' : 


COS 0' 

77 


f 


cos u ' cos iu ' du ' 

( 1 - sin 0' cos u ' ) 2 


i tan 1 1 


1 , 

2 ^ 


2 cos 2 


1 

2 0' cos 0' 


tan 


+ 


2 0* sin 0' 
cos 2 0 1 


(143b) 


But, multiplying Equation (143a) by: 


1 - sin 0' cos u ' 
1 - sm 0* cos u ; 


gives: 


77 


f 

Jo 


cos iu ' du ' 

( 1 - sin 0' cos u ' ) 2 


sin 0' 

77 


f 


cos u ' cos i u ' du ' 

( 1 ~ sin 0' cos u ' ) 2 


tan 1 2 0 # 
COS 0' 


Equations (143b) and (143c) give: 

i 


I 

77 


1' 

Jo 


cos iu * du ' 

( 1 - s in 0' cos u ' ) 2 


I 

I 

Applying this procedure a second time: 


tan 1 2 0 ' 
COS 3 0' 


( 1 + i cos 0' ) 


77 


f 


cos i u ' du ' 

( 1 ~ sin 0' cosu') 3 


tan 1 “2 0' 

“ ]2 + sin 2 0' + 3i cos 0' + 

2 cos 5 0' L 



(143c) 


(143d) 


(144) 


Applying the above results to Equations (141), we have: 



tan 1 1 "" 2 " 0 * 


|sin 2 0' + i cos 0' + i 2 cos 2 0'j 

(145a) 

Mi' 

9 1 , - 
4 cos 2 0 cos- 3 



■ _i 1 , 

tan '2 0 


i + i 2 co s 0 ' J 

(145b) 

v 

i 

2 1 , , 
4 cos z 2 " 0 cos 

0' 


Consequently, the perturbation equations for the second part of the perturbation function, may 
be written as: 

~gjjr - r 3 P/i 0 + r 3 cos u ' + r 3 P/7 2 cos 2u' + *** 

+ r 3 Qj sin u' + r 3 Q i^ 2 sin 2u' + (146a) 
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dk 


r 3 PjM 0 + r 3 Pj cos u ' + r 3 P t /i 2 cos 2u ' + *■* 



+ r 3 Qj i'j sinu' 

+ r 3 Qj v 2 sin 2u ' 

+ 


(146b) 

dE 






dif _ 

r 3 P 2 A^o + r 3 P 2 ^ 1 cosu' 

+ r 3 P 2 (jl 2 cos 2 u ' 

+ 




+ r 3 Q 2 v x sin u ' 

+ r 3 Q 2 v 2 sin 2 u ' 

+ 


(146c) 

dp _ 
dk 

r 3 P ^ Mq + r 3 P 3^1 cosu' + 

r 3 P 3 M 2 cos 2u ' + • 

-n 






cos i 

(146d) 


+ r 3 Q 3 v j sinu ' + 

r 3 Q 3 v 2 sin 2u ' f- • 

- 



Art 

r 3 P 4 fj . 0 + r 3 P 4 /Xj cos u' + 

r 3 P 4 fji 2 cos 2u ' + • 

-i 



nq 

dk 




cos i 

(146e) 


+ r 3 Q 4 ^ j sinu' + 

r 3 Q 4 v 2 sin 2u ' + • 





where: 


except: 


m ' 



m' sin 0' 
a ' 3 cos 3 0' 


(147) 


The quantities P, Q, P lf Q 1? etc., are given by Equations (139). The total value of dy/dk, d^/dk, 
etc., is given by the sum of Equations (138) and (146). 

The series given by Equations ( 146 ) may not converge when r > a unless only a very small 
segment of the ellipse near aphelion is represented by the superior anomaly. This results since 
r 3 is a factor and if it is large, the convergence is weak at best. To overcome this difficulty, we 
divide the perturbation function as follows: 



m ' 

"l r ' 


\ 

r ' r 

n 

“A " m ' 

f + 7» H J 

+ m ' 

r 

1 H - — H 
r 1 r 


(148) 
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The first term has already been considered. The third term may be integrated in finite terms and 
will be considered later. 

Consider the second term: 


Then: 



dn d H 

d f ' m ' r 2 di 


, ill iH 

“ m r 3 df 


(148a) 


dO r' m' 

r dr D 2m 2 H+ r 


rr r * 

2m - H + m' r 
r J r 3 


(148b) 


en 

dH 


(148c) 


Equations (148) will be identical with Equations (131) if, in Equations (131), A 3 is replaced by 
“r 3 , and -3m' rr'/r 2 is added to the second Equation (131). Proceeding as we did in deriving 
Equations (138), we find: 


dy 

dk 


m ' R + e ' 


6 m' a ' a 2 y 
r 


ndt 
r 2 dk 


sin f sin ( f + T) 


- 2m' P- 


6 m ' a ' a 2 y 
r Vl - e 2 


ndt 
r 2 dk 


sin f sin ( f + T) 


6 m' a' a 2 y' rl- e ' 2 ndt 
2- Q ^ Sinf ,in(f + n 


(149a) 


d! 

dk 


. . 6 m' a' a 2 7 ndt „ . . „ 

m R, - e' ^ 2 * 2 ~ cos f sin (f + P) 


il 


r 2 dk 


6 m' a' a 2 7 ndt 

2m' P. + r , •• — T~ cos f sin ( f + T) 

1 r i\ - e 2 r 2 dk V J 


2m' Qj 


6m 1 


a 2 7 ‘ 






ndt 
r 2 dk 


co s f sin ( f + T' ) 


(149b) 
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ds 

dk 


dp 

dk 


dg 

dk 


R 2 + 2m' P 2 cos u' 

+ 2m ' Q 2 sin u 'J 

(149c) 

R 3 + 2m' P 3 cos u' 

+ 2m' Q 3 sin u 'j cos i 

(149d) 

R 4 + 2m' P 4 cos u ' 

+ 2m' Q 4 sin u 'J cos i 

(149e) 


Equations (149) will be very convergent when the superior anomaly is used but will not be 
convergent with the inferior anomaly. Consequently, Equations (146) are used with the inferior 
anomaly and Equations (149) with the superior anomaly. 

Consider the third term of Equation (148) 


.0 


3 



Proceeding as in Equations (131), (132), (133), (134), (135), (136), and (137), the perturbation equa- 
tions for n 3 are: 


A 1 

dt 


2m' 


an 


cos f cos ( f + P) cos 


f ' 


cos f + e 
1 - e 2 


cos ( f + T) cos f ' 


r r 


ar a 

- 2y — r sin f sin ( f + P) cos f ' - y cos P cos f ' 
r J v 1 


r r 


, cos f + e „ , . - 

+ 7 2 — cos (f + r ) sin f ~ - “ 

1 - V r 1 r 


+ y ' — “ cos f cos(f + P') sin f' 


- 7* — — cos P 1 sin f ' 


ar ' a 

” 27 ' „ sinfsin(f + P')sinf' " , sinf 

r r * 


(150a) 


52 



dt 


2m' 


an 


ar ' 



sin f cos ( f + P) cos f ' 


sin f _ 

+ 7 : — * cos ( f + r> cos f ' 
1 “ e 



, sin f 

+ 7 j _ 2 cos ( ^ + T' ) sin f ' 





sin f cos (f + r') sin f ' + y' sinP' sin f ' 


+ 2y' j cos f sin ( f + P' ) sin f ' 



cos f 


(150b) 


dg 

dt 



y cos ( f + P) cos f ' 




+ y' cos (f + T') sinf' 




(150c) 


_d£ 

dt 


an 

il - e 2 


co s i sin I 


y' sin G ' sin f cos f ' 




“ y sin G sin f sin f ' 




(150d) 


dq 

dt 


/ an 

m n « cos i sin I 

Tl - e 2 


y ' sin G ' cos f cos f ' 



- y sin G cos f sin f ' 




(150e) 


The perturbation equations have now been written in a form proper for integration in all cases. 
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V. SERIES EXPANSION OF A" 


Consider a general periodic function. 

X - a 0 + a i cos x + a 2 cos 2 x + ■ • • + a n cos n x 

+ sin x (/3 0 + ySj cos x + /3 2 cos 2 x + • • • + /S n _j cos n 1 x) (151) 

We will show that this can always be written in the form: 

X - c [l - qj cos (x - Qj)J |l - q 2 cos (x - Q 2 )J ••• [l - q n cos (x - Q n )J (152) 

Define: 

tan ~2 x = y (153) 


Then: 


2y 
1 + y s 


l-y 2 

1 + y 2 


COS X 


(154) 


Equation (151) may be written as: 


(l + y 2 ) n 


a 0 ( 1 +y 2 ) n + a, (l-y 2 ) (l+y 2 )"" 1 


a 2 (i ~ y 2 ) 2 (l + y 2 ) n 2 + ■ • • + a n (l _ y 2 ) n 


2y 


(1 +y 2 ) n 


0 O (i +y 2 ) n ‘ 1 + 0 t (l-y 2 ) (l +y 2 )"' 


0 n -i (i-y 2 )"' 1 


Suppose: 


X(l+y 2 )” = 0 

Equation (156) is of degree 2n and may be written: 


(155) 


(156) 


0 = C 2n y 2 " + C 2n-I y2n "' + ^ y t C Q 


(157) 
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Let the roots be: 


Then: 


3 1- a 2> ••• • a n 


b l- b 2’ • b n 


X(l+y 2 ) n = c 2n (y-a,) (y-bj (y-a 2 )(y-b 2 ) ••• (y-aj (y-b n ) 


(158) 


or: 


where: 


X(l+y 2 )" = C 2n (y 2 -s iy + Pl ) (y 2 -s 2 y + p 2 ) ( y 2-s n y + pJ 


a l + b l = s l 


a l b l - P, 


a, + b 2 = s 2 a 2 b 2 = P 2 


a +b - s ab = p 

n n n n n r n 


We can always select s ; and p, to be real. Since: 


(159) 


(160) 


1 + y 2 


tan -y x 


2 i 

sec 4 x 


Consequently: 


,2n 


1 / 2 1 1 \ 

2 x I tan 4 ~2 x - s x t an y x + P x ] 

/ 2 1 1 \ / 2 1 1 
X It an 4 2 x - s 2 tan yx + pJ * * * tan^ 2 x- s n tan ^ x + P n 


(161) 


If we multiply each factor by cos 2 1/2 x, it will take the form: 


2 I I i. ,1 

sin 2 x “ s sin 2 x cos 2 x + P cos *2 x 
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or: 


Define: 


and: 


y (i +p) - y(i 

~ p) cos x 

1 

- y s sinx 

1 ~ p i 

l+p, s q,cosQ, 

s l 

1 + Pj 

QjSinQj 

i -p 2 

i+p 2 - q 2 cos Q2 

__ S 2_ 

1+P 2 

= q 2 sin Q, 

1 - Pn 

S 

ii 


i+ Pn = <i„ cos Q„ 

1 + Pn 

q sin Q 

Ml x 1 


C 5 ^ C 2n (1 + Pl) (* + P 2 ) ••• (* + Pn) 


or 


C 0 l 1 + Pt) l 1 + P ; ) l 1 + P 3 ) ' • ' I 1 + Pn) 

2 " a l a 2 a 3 a n b l b 2 b 3 ••• b n 


Equation (155) becomes: 


X 


[l - q, cos (x - Qj)] 



q n c°s 


i* 



that is, we have proved Equation (152). 

The case n = 2 will of some interest. Comparing terms in Equations (155) and (157), 


(a 0 - a] + a 2 ) y 4 + 2(/3 0 - /3i) y 3 + 2(a 0 - a 2 ) y 2 + 2(/3 0 + /3,) y + (a 0 + a, + a 2 ) - 0- 

Equation (151) is: 


X ' a 0 -a 1 cosx + a 2 cos 2 x-/3 0 sinx + /3 1 sinxcosx 


(162) 


(163) 


we have: 
(164) 


(165) 
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Write: 


X 


1-^ 


q . 

~ cos(x 





Comparing terms, we have: 


qq! 

a o “ c + sin Q sin Q 2 

a i = q cos Q + q 1 cos Qj 
QQi 

a 2 = ~ cos (Q + 0i) 

/3 0 = q sinQ + qj sinQ; 

qq i / 

~ ~ sin^Q + Qj) 

Hence: 

, P l 

tan (Q +< ?i) = ^7 

Define: 

0 + Q 1 - 20 

Q - Qj s 

Then: 

tan 2^ - — 

tt 2 

Q = 0 + oo 

Qj ° <9 - o) 


(166) 


(167) 


(168) 


(169) 
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Eliminating Q and Qj from Equations (166), we have: 


Define: 


c + „ (sin 


2 d~ 


a 2 - ~ — cos 2d 




^q + q^ cos d cos * - (q-q^ sin d sin .«.• 


(q + Qj) sin 0 cos i, + - qj) cos d sin < 


qqi 

sin 2d 


Uq a Q + a 2 sin 2 fc> - j3 j sin 0 cos 6 

a j' - Uj cos 0 + fi Q sin U 
a 2 ' a 2 cos + (i x sin 2lf 

/S 0 ' = /3 0 cos - a j sinf^ 


(170) 


(171) 


The quantity ^ is determined by Equation (169). With the aid of Equations (171), Equations (170) 
become: 


a . ^ - (q + q ^ cos u: 

, - 

a 2 C 

/V 7 (q-qi) sin,< ' ( 172 ) 

We are going to require that the q f s always be positive. In this particular case, this requires 
that ■ be positive. Consequently, we select the value of o satisfying Equation (169) so thata 2 ' is 



positive. Define: 


q cos 01 - cl i - r\ qjcosoj - rj 

qsino; ^ / 3 Q ‘ + g q i sin ox = £ 

qq, 

c s a; + 5 or i - sin 2 co (173) 

Equations (172) and (173) give: 

K + 0 S = (V + ?) € (174a) 

“a (“o' + o = (»,' - v) v + (/y + i) s (174b) 

[ a i'-o)€ = (V + ?) r i (174c) 

Equations (174) must be solved to obtain i, v, and ?. Equations (174a) and (174b) give, by 
subtraction: 


K' + 0 K' " 0 D 


Equation (174a) may be written: 


4 @0 2 + ( a o* + £) £ 


= (y/v + ?) 2 


Using Equation (174), we have: 


2 4/ + 


- i ,1 

' 2 “l T) 


Equations (176) and (177) yield: 


4 + ( a o + 0 £ r 4 a o 


-i„ .2 if 


Now, Equations (174a) and (175) may be written as: 


(“o' + 0 £ 

A> + £ ~ | " 


a , - V 


(“o' + 0 K "£) 


(175) 


(176) 


(177) 


(178) 
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Inserting these into Equation (174c) gives: 


(179) 


Eliminating : 2 /V between Equations (178) and (179) gives: 

:* ' • (l ■/* •iv’--.’-.) - 


(180) 


Equation (180) determines the possible values of || Equations (174a) and (175) may be written in 
the following form. 


- ( + ■' ) l - 0 

(181) 

o' + 4) (" 2 ' - 0 

(182) 


For each value of we can determine values of ~ and tj. To determine which roots of Equations 
(181) and (182) are to be taken together, the following condition must be satisfied. 


« i ~ 7 > IV 4 -f 
n •: 


(183) 


Which pair of roots satisfying Equation (183) is taken is unimportant as this simply results in a 
permutation of the factors in X. 

We may have some use of the particular case, B l = 0, and i 2 very small. If = 0, Equa- 
tion (169) gives o = 0. Equations (171) become: 


a _ - a „ 


$ o ^0 


Also, from Equation (169): 


Define: 


» 77 0 7 ~ 0, 


f sin F 
n j f cos F 
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Equations (180), (181), and (182) become: 


4a 


£ - a 2 sin 2 F + 


0 a 2 4 (a 0 a 2) 


r - — 1 — l 2 - — r 3 

f 2 ^ f2 S f2 ^ 


f - ( a ° + L 1 4 * 2 
^ " /So ^ 


(“o + t)k-0 1 


+ T7 2 


If a 2 is small, these equations may be quickly solved by iteration. Let us now return to the gen- 
eral case. 

Recall Equation (155). 

X = c [l - q, cos (x - Q,)] ••• [l - q n cos (x - Q n )] 


Let: 


f. - c - q. cos0. 

<t>i =-■ x - Q. 


In problems of interest to us, q. < 1 . Define: 


<t> l = 180° - 20/ 


! Then: 

j 

i 

I 




f . - 

i 


(c + qj 1 


c + q . 


sin <p 


i 


(184) 


(185) 


I We may write: 

► 

► 

r 


X " c • f, • f 2 ••• f n 


We will be required to determine X’ n/2 where n is an integer. Now: 


V~ n/ 2 m “ n/2 c - n / 2 . f - n/2 ... f - n/ 2 

A C r i r 2 r n 
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But, f . is given by Equation (185) and this is of the form of Equation (27a) which we have developed 
in detail for n = 1. In the general case, the equations of section n become: 

(l - e 2 sin 2 k) _n/2 = a 0 - 2a 2 cos 2k + • • - 


a 


2 i 



cos 2ik dk 
(l - e 2 sin 2 kp 2 


0 - (2i+n-2)e 2 a 2 ._ 2 -4i(2-e 2 )a 2i +(2i-n + 2)e 2 a 2i + 2 


2i + n ~ 2 6 2 _ a 2 i 

4i 2 -e 2 721 S a 2 i -2 

(2i -n) (2 i + n - 2) / e 2 V 
16 i ( i ~ 1) \2 ~ e 2 ) 


The approximation for y 2i is: 


7 


2 i 


: 2 \ X 


n (n - 2) i_ 

4i( i + 1) sin 2 X 


where: 


e - sin 0 
cos - tan 

We will be somewhat interested in the case n = 3. We have: 


(45 - \ x) 


a < 3 > 
a o 


-*r 


dk 

(l - € 2 sin 2 k) 3/2 


Since: 


(i) 


and: 


-If 


dk 


e 2 sin 2 k) 


1/2 


a < 3 > 

u 2 77 


d 


cos 2k dk 


0 (l-e 2 sin 2 k) 3/2 
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we have: 


a <*> 

0 


(l - | c 2 ) 


a 0 (3) 


+ 2 


e 2 a 2 < 3 ) 


But, the definition of p 2 . gives: 


( 3 ) = 


P 2 <3) a 0 <*) 


and 


*o (3) 


a 0 ^ 


- 2 e 2 + 7 e 2 p 2 (3) 


Returning to the expansion of f n/2 , if we define: 


sin x = ■£- 


(186) 


we write: 


f-n/2 


M 0 + 2/u, cos 4> + 2/J.J cos 20 + • • • 


(187) 


where: 


‘ ^0 * P2 * P4 P 2 i 


(188) 


For n = 1; 




(i) 


( c + <3i) 


“1/2 rn 


f 

Jo 


dk 

2 q A \ 1/2 

1 ■ 7T^7 sin 2 0.'l 


(189) 


For n = 3: 


( 3 ) 


( C + ^i) 


- 3/2 


dk 


2q \ 3/2 

1 " c + q. si " 2 <t>i 


( 190 ) 
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Hence: 


(i) 


( 3 ) 


M o 
c - q p 


( 3 ) 


We now write: 


|l - q J cos (x - Qj) J n/2 - M 0 (1) + !) cos Qj cos x + 2/x 2 (1) cos 2Qj cos 2x 


2/x/ 1 * sinQj sin x * 2^ 2 (1) sin 2Qj sin 2x 


- q 2 cos (x - Q 2 )J n/2 - m 0 (2) + 2/Xj (2) cos Q 2 cos x + 2/x 2 (2) cos 2Q 2 cos 2x 


+ 2^j (2) sin Q 2 sin x + 2^ 2 (2) sin 2Q 2 sin 2x + • • ■ etc. 


write: 


* n/2 


[m, 0] + 2[m, 1] cos x + 2Lm, 2] cos 2x + 


+ 2[m, 1] sin x + 2[m, 2] sin 2x f 


where, m is the number of factors in X. Then: 


[m, l\ 




-n/2 > „(1) . .,(2) ... 


^iV 


Wm) cos(i' Q 1 + i"Q 2 4 ... + i (m) Q m ) 


and: 


lm, l] 


-n/2 


^~ l m/, 0 • m/ 2) ••• sin(i' Q,+ i"Q 2 + ■•■<■ i <m ^ 0 m ) 


where: 


l = i ' + i" + • • * t i < ni > 


The sum is over all values of combinations of i ' , i ", • • • , i (m) , which sum to 


( 191 ) 


.( 192 ) 


( 193 ) 


( 194 ) 


( 195 ) 


64 



Recall the elliptic equations. 


r ' - a ' ( 1 - e ' cos u ' ) 


r ' cos f ' - a' (cosu' - e' ) 


r' sinf' n a' Yl - e ' 2 sinu' 


Equations (126) and (127) were: 

H = cos ( f + v + k) cos ( f ' + v - k) + cosIsin(f+^+k)sin(f'+v-k) 


and: 


A 2 = r 2 + r ' 2 - 2rr ' H 

Combining these with Equations (134) for y, G, r, and y', G' , r', we have: 

r ' H ' Aa ' cos u ' + Ba ' sin u ' - Ae ' a' (196) 


where: 


A = y sin ( f + f) 


B = y‘ Yl - e'~^ sin ( f +T') (197) 


and: 


A 2 - a' 2 + r 2 + 2e'Aa'r - 2 (Aa ' r + a' 2 e') cosu' - 2Ba ' r sinu' + a' 2 e' 2 cos 2 u' (198) 

We have now expressed A 2 as a periodic function of u'. By means of the equations of elliptic 
motion, we may convert this to an infinite series in powers of sines and cosines of the mean anomaly. 
If, in this series, we neglect powers higher than the nth, we will have A 2 as a periodic function of 
sin g ; and cos g' of the nth degree and the results of this section may be applied. The actual appli- 
cation will be discussed in Section VI. 
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VI. A DISCUSSION OF THE APPLICATION 
OF THE METHOD OF PARTIAL ANOMALIES 


Sections I- V of this report provide most of the theory of the method of partial anomalies neces- 
sary to be able to apply it to problems of interest. This section will summarize the earlier results 
and fill in the few gaps left in the previous sections. 

The method of partial anomalies as developed by Hansen (Reference 1) applies directly to the 
perturbations produced by a body, considered as a point mass, on the unperturbed two body elliptic 
motion. The generalization to more than one perturbing body is obvious, since the perturbations 
may be considered as additive in all cases likely to occur. An important application of this method 
may be to high eccentricity artificial Earth satellites as well as to interplanetary probes. In fact, 
this method will converge for any value of e, the eccentricity. The only difficulty occurs if the 
perturbation is so large that the orbital character of the motion is completely changed so as to 
reverse the direction of motion. Needless to say, this causes difficulty in any known method. 

Consequently, in order to apply the method to artificial satellites, the method must be adapted 
to the perturbing force due to the departure of the Earth r s shape from sphericity. The author of 
this report is currently beginning an attack on this problem in order to determine the quantitative 
value of this method as applied to high eccentricity satellites and probes. Let us now return to 
the problem at hand. 

Hansen was interested in applying this method to comets, most of which have high eccentrici- 
ties, and partially determined the perturbations on Encke’s comet produced by the Earth. In what 
follows, I shall refer to the comet and the Earth, rather than the perturbed body and the perturbing 
body, respectively. It is clear that this ’’specialization” does not cause a loss of generality while 
it may maintain Hansen’s line of thought in a clearer fashion. 

The problem may be divided into three major parts. First, the definition of the inferior and 
superior anomalies and their auxiliaries, K 1? K 2 , k 3 , K 4 , and the intermediate anomalies. The 
equations of elliptic motion are then expressed in terms of the appropriate partial anomaly. This 
part of the problem is dealt with in Section I. We have seen that the comet’s ellipse may be divided 
into any number of segments in an arbitrary manner. Exactly how this is done will be determined 
by the specific problem and by the choice between rapid convergence and the number of series 
which must be developed. Since a new set of series is required for each segment of the ellipse, 
it is obvious that we will reach the point where the work required by n sets of series is not worth 
the greater convergence obtained by dividing the ellipse into n segments rather than n - 1 segments. 

Once we have expressed the motion of the comet in terms of the partial anomalies, we turn to 
the second part of the problem. This is the development of various functions, such as the perturb- 
ing function, in trigonometric series of multiples of sines and cosines of the partial anomalies. 

Once expressed in this form, integrations are readily performed. This development in trigonometric 
series has been discussed in Sections II, IV, and V. 
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The expressions for sin f and ndt in terms of the partial anomalies in Section I contain radicals 
such as (l - e 2 sin 2 k)~ 1/2 In Section II, these radicals are developed into trigonometric series con- 
venient for integration. The end result of Section II is that the elliptic motion of the comet has been 
expressed in terms of the partial anomalies in a readily integrable form. 

In Section IV, the perturbation equations are developed into trigonometric series involving the 
eccentric anomaly of the Earth. The sines and cosines of this eccentric anomaly may be converted 
to sines and cosines of the mean anomaly of the Earth by means of the elliptic equations and is 
generally performed by means of Bessel functions. This is discussed in standard texts on celestial 
mechanics. Finally, by means of Equation (84): 

g ' = n ' t + c ' 

we express the perturbation equations in terms of sines and cosines of c'. Since c' is independent 
of the partial anomalies, integrations are readily performed. 

In Section V, A' n is developed into trigonometric series involving the Earth's eccentric anomaly. 
Again, these may be converted to a trigonometric series in c' and again we have a form convenient 
for integration. 

The second part of the problem having been completed, we turn to the third and final problem. 
This is the determination of the constants of integration, which was discussed in Section III. Since 
each segment of the ellipse is represented by a different set of series, each segment will have its 
own constants of integration. These are determined by initial conditions, such as the initial values 
of the orbital elements of the comet, and by the condition that we have continuity at the points of 
separation which divide the various segments of the ellipse. 
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